Representation, searching and discovery of patterns of bases in complex RNA structures

We describe a graph theoretic method designed to perform efficient searches for substructural patterns in nucleic acid structural coordinate databases using a simplified vectorial representation. Two vectors represent each nucleic acid base and the relative positions of bases with respect to one another are described in terms of distances between the defined start and end points of the vectors on each base. These points comprise the nodes and the distances the edges of a graph, and a pattern search can then be performed using a subgraph isomorphism algorithm. The minimal representation was designed to facilitate searches for complex patterns but was first tested on simple, well-characterised arrangements of bases such as base pairs and GNRA-tetraloop receptor interactions. The method performed very well for these interaction types. A survey of side-by-side base interactions, of which the adenosine platform is the best known example, also locates examples of similar base rearrangements that we consider to be important in structural regulation. A number of examples were found, with GU platforms being particularly prevalent. A GC platform in the RNA of the Thermus thermophilus small ribosomal subunit is in an analogous position to an adenosine platform in other species. An unusual GG platform is also observed close to one of the substrate binding sites in Haloarcula marismortui large ribosomal subunit RNA

A long-acting GH receptor antagonist through fusion to GH binding protein.

Acromegaly is a human disease of growth hormone (GH) excess with considerable morbidity and increased mortality. Somatostatin analogues are first line medical treatment but the disease remains uncontrolled in up to 40% of patients. GH receptor (GHR) antagonist therapy is more effective but requires frequent high-dose injections. We have developed an alternative technology for generating a long acting potent GHR antagonist through translational fusion of a mutated GH linked to GH binding protein and tested three candidate molecules. All molecules had the amino acid change (G120R), creating a competitive GHR antagonist and we tested the hypothesis that an amino acid change in the GH binding domain (W104A) would increase biological activity. All were antagonists in bioassays. In rats all antagonists had terminal half-lives >20 hours. After subcutaneous administration in rabbits one variant displayed a terminal half-life of 40.5 hours. A single subcutaneous injection of the same variant in rabbits resulted in a 14% fall in IGF-I over 7 days. IN CONCLUSION: we provide proof of concept that a fusion of GHR antagonist to its binding protein generates a long acting GHR antagonist and we confirmed that introducing the W104A amino acid change in the GH binding domain enhances antagonist activity

COGNAC: a web server for searching and annotating hydrogen-bonded base interactions in RNA three-dimensional structures

Hydrogen bonds are crucial factors that stabilize a complex ribonucleic acid (RNA) molecule's three-dimensional (3D) structure. Minute conformational changes can result in variations in the hydrogen bond interactions in a particular structure. Furthermore, networks of hydrogen bonds, especially those found in tight clusters, may be important elements in structure stabilization or function and can therefore be regarded as potential tertiary motifs. In this paper, we describe a graph theoretical algorithm implemented as a web server that is able to search for unbroken networks of hydrogen-bonded base interactions and thus provide an accounting of such interactions in RNA 3D structures. This server, COGNAC (COnnection tables Graphs for Nucleic ACids), is also able to compare the hydrogen bond networks between two structures and from such annotations enable the mapping of atomic level differences that may have resulted from conformational changes due to mutations or binding events. The COGNAC server can be accessed at http://mfrlab.org/grafss/cognac

Evidence that a salt bridge in the light chain contributes to the physical stability difference between heavy and light human ferritins.

Human ferritin, a multimeric iron storage protein, is composed by various proportions of two subunit types: the H- and L-chains. The biological functions of these two genic products have not been clarified, although differences in reactivity with iron have been shown. Starting from the hypothesis that the high stability typical of ferritin is an important property which may be relevant for its iron storage function, we studied ferritin homopolymers of H- and L-chains in different denaturing conditions. In addition we analyzed 13 H-chain variants with alterations in regions conserved within mammalian H-chains. In all the denaturation experiments H-chain ferritin showed lower stability than L-chain ferritin. The difference was greater in guanidine HCl denaturation experiments, where the end products are fully unfolded peptides, than in acidic denaturation experiments, where the end products are peptides with properties analogous to "molten globule." The study on H-chain variants showed: (i) ferritin stability was not affected by alterations of regions exposed to the inner or outer surface of the shell and not involved in intra- or inter-chain interactions; (ii) stability was reduced by alterations of sequences involved in inter-subunit interactions such as the deletion of the N-terminal extension or substitutions along the hydrophobic and hydrophilic channels; (iii) stability was increased by the substitution of 2 amino acids inside the four-helix bundle with those of the homologous L-chain. One of the residues is involved in a salt bridge in the L-chain, and we concluded that the stability difference between H- and L-ferritins is to a large extent due to the stabilizing effect of this salt bridge on the L-subunit fold

Computation of protein geometry and its applications: Packing and function prediction

This chapter discusses geometric models of biomolecules and geometric constructs, including the union of ball model, the weigthed Voronoi diagram, the weighted Delaunay triangulation, and the alpha shapes. These geometric constructs enable fast and analytical computaton of shapes of biomoleculres (including features such as voids and pockets) and metric properties (such as area and volume). The algorithms of Delaunay triangulation, computation of voids and pockets, as well volume/area computation are also described. In addition, applications in packing analysis of protein structures and protein function prediction are also discussed.Comment: 32 pages, 9 figure

A Burkholderia pseudomallei toxin inhibits helicase activity of translation factor eIF4A

The structure of BPSL1549, a protein of unknown function from Burkholderia pseudomallei, reveals a similarity to Escherichia coli cytotoxic necrotizing factor 1. We found that BPSL1549 acted as a potent cytotoxin against eukaryotic cells and was lethal when administered to mice. Expression levels of bpsl1549 correlate with conditions expected to promote or suppress pathogenicity. BPSL1549 promotes deamidation of glutamine-339 of the translation initiation factor eIF4A, abolishing its helicase activity and inhibiting translation. We propose to name BPSL1549 Burkholderia lethal factor 1

X-ray analysis of protein dynamics

The structure of human lysozyme has been refined at 1.5 Å resolutions in terms of 1026 protein atoms and 130 water molecules. The current R-factor is 0.187 for the 19000 observed reflections. Although analysis of the structure of the enzyme yields important results, analysis of the vibration data seems even more important. A detailed comparison with the results obtained from and independent refinement of hen lysozyme at 2 Å resolution shows that both molecules have a very similar pattern of atomic displacements which clearly reflects the molecular structure of the protein. The region of high motion includes part of the active site and those residues that undergo conformational changes on substrate binding. These results suggest that analysis of protein dynamics is possible with X-ray diffraction studies

Selection of reagents for combinatorial synthesis using clique detection

The identification of a database subset in which all pairs of molecules have similarities less than some user-defined threshold is shown to be equivalent to the graph-theoretic problem of clique-detection. A comparison of several different clique-detection algorithms demonstrates that the algorithm due to Babel is by far the most efficient for this particular problem. The database subsets identified by the Babel algorithm are slightly more diverse than those identified by alternative, less time-consuming algorithms for dissimilarity-based compound selection

Crystal structure of DNA recombination protein RuvA and a model for its binding to the Holliday junction.

The Escherichia coli DNA binding protein RuvA acts in concert with the helicase RuvB to drive branch migration of Holliday intermediates during recombination and DNA repair. The atomic structure of RuvA was determined at a resolution of 1.9 angstroms. Four monomers of RuvA are related by fourfold symmetry in a manner reminiscent of a four-petaled flower. The four DNA duplex arms of a Holliday junction can be modeled in a square planar configuration and docked into grooves on the concave surface of the protein around a central pin that may facilitate strand separation during the migration reaction. The model presented reveals how a RuvAB-junction complex may also accommodate the resolvase RuvC

Structural similarities between Escherichia coli RuvA protein and other DNA-binding proteins and a mutational analysis of its binding to the Holliday junction.

Comparison of the structure of Escherichia coli RuvA with other proteins in the Protein Data Bank gives insights into the probable modes of association of RuvA with the Holliday junction during homologous recombination. All three domains of the RuvA protein possess striking structural similarities to other DNA-binding proteins. Additionally, the second domain of RuvA contains two copies of the helix-hairpin-helix (HhH) structural motif, which has been implicated in non-sequence-specific DNA binding. The two copies of the motif are related by approximate 2-fold symmetry and may form a bidentate DNA-binding module. The results described provide support for the organization of the arms of the DNA in our RuvA/Holliday junction complex model and support the involvement of the HhH motifs in DNA binding